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I consider a Josephson junction crossing the middle of a thin rectangular superconducting strip 
of length L and width W subjected to a perpendicular magnetic induction B. I calculate the 
spatial dependence of the gauge-invariant phase difference across the junction and the resulting B 
dependence of the critical current Ic{B). 
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I. INTRODUCTION 

Grain-boundary Josephson junctions play an impor- 
tant role in thin films of YBasCuaOr-a (YBCO).!"'' Var- 
ious theoretical approaches have been taken to under- 
stand the physics of Josephson junctions in thin films. 
When the film thickness d is less than the London pene- 
tration depth A, the current density j is practically uni- 
form across the thickness, and the characteristic length 
governing the spatial distribution of the magnetic field 
distribution is the Pearl length,® 

A = 2AVd. (1) 

However, various studies have shown that there is a non- 
local relationship between the Josephson-current distri- 
bution in the vicinity of a Josephson vortex core and the 
magnetic field these currents generate, ^"^^ and the char- 
acteristic length describing the spatial variation of the 
gauge-invariant phase across the junction is (in SI units) 

t = 0o/4^MoA2ic, (2) 

where <j)Q = h/2e is the superconducting flux quantum 
and jc (assumed to be independent of position) is the 
maximum Josephson current density that can flow as a 
supercurrent through the junction. 

The integral equations relating the gauge-invariant 
phase difference across the junction to the magnetic field 
generated by a vortex in a long Josephson junction in a 
thin (d < A) film of lateral dimensions large by compari- 
son with A and i were examined analytically and solved 
numerically in Ref. 6 for arbitrary ratios oii/ K. The case 
of a short Josephson junction bisecting a long supercon- 
ducting strip of width W was studied in Ref. 7 under the 
assumptions that W <^ K and W ^ £. 

In this paper I revisit the latter problem by consider- 
ing a thin rectangular uniform superconducting strip of 
length L, width W, and thickness d {d < X) divided into 
two halves by a Josephson junction at a; = 0, as shown in 
Fig. 1. An applied magnetic induction B — zB induces 
screening currents in the film. However, I consider here 
only the simplest case for which A = 2A^/d is much larger 
than the smaller of L and W , such that the self-field gen- 
erated by the screening currents can be neglected.'' The 
purpose of this paper is to calculate how the screening 
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FIG. 1: Sketch of the thin rectangular superconducting strip 
of width W and length L considered here. An in-line Joseph- 
son junction of width d; (bold line) is at the center at x = 0. 
A magnetic induction B is applied in the z direction. Current 
leads of separation L (not shown) symmetrically feed current 
/ to the sample along the x direction. 

currents induced in response to B affect the B depen- 
dence of the maximum Josephson critical current Ic{B). 

In Sec. II I give a brief discussion of the derivation of 
the basic equation for the gauge-invariant phase differ- 
ence A7(y) across the junction, in Sec. HI I present the 
solutions for A7(y) for arbitrary ratios of L/W , and in 
Sec. IV I briefly summarize the results. 



II. GAUGE-INVARIANT PHASE DIFFERENCE 

In the context of the Ginzburg-Landau (GL) 
theory, ^^'^^ the superconducting order parameter can be 
expressed as ^ = where 5*0 is the magnitude 

of the order parameter in equilibrium, / = |\['|/5'o is the 
reduced order parameter, and 7 is the phase. The second 
GL equation (in SI units) is 

j--^(^ + |^V7), (3) 

where A is the vector potential and B ~ V x A \s the 
magnetic induction. Since j is a gauge-invariant quan- 
tity, so is the quantity within the parentheses on the 
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right-hand side. Different choices for the gauge of the 
vector potential A resuh in different expressions for 7. 

With a sinusoidal current-phase relation, the Joseph- 
son current density in the x direction across the junction 
of width dj at X = is jx{y) = jc sin Aj{y), where jc 
is the maximum Josephson current density and Aj{y) is 
the gauge-invariant phase difference between the left (a) 
and right (b) superconductors, 



90 J-di/2 



A^{x,y)dx. 

(4) 

I assume here that the induced or applied current den- 
sities ja and jh on the left- and right-hand sides of the 
junction are so weak that the suppression of the mag- 
nitude of the superconducting order parameter is negli- 
gible, such that / = 1. A simple relation between these 
current densities and the gauge- invariant phase difference 
can be obtained by integrating the vector potential A 
around a very narrow rectangular loop of width di in the 
xy plane that just encloses the junction (with the bottom 
end at the origin and the top end at y), neglecting the 
magnetic flux up through the contour, making use of Eq. 
(3) with / = 1 for those portions of the integration along 
the sides of the junction, and noting that, by symmetry, 
ia!/(0,y) = -jby(0,y): 



Aj{y) = A70 + 



jby{0,y')dy\ (5) 



where A70 = A7(0), such that 



dAj{y)/dy = {AwiioX^ /(j)o)jby{0,y)- 



III. SOLVING FOR A7 
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I next assume that the Josephson coupling is so weak 
that the currents ja and jb on the left- and right-hand 
sides of the junction induced in response to the applied 
magnetic induction B are far larger than the Josephson 
current density. This is equivalent to the assumption that 
W <^ £. Since easily can be obtained by symmetry 
from jb, I calculate only in the region a; > and 
suppress the subscript b. 

With the gauge choice A = —xBy, since V • j = [see 
Eq. (3)], V^7 = must be solved subject to the boundary 
conditions following from jx{0,y) = jx{L/2,y) = and 
jy{x,±W/2) = 0, 

7.(0, y) = y) = 27rBy/<^o, (7) 

7,(x,±W^/2) = 0, (8) 

where 7^ = dj/dx and 7^ = d^/dy. The solution, ob- 
tained by the method of separation of variables, is (up to 
a constant) 
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FIG. 2: Contour plot of the stream function S{x,y) in the 
right half {x > 0) of a strip of length L = 2W. The contours 
correspond to the streamlines of the current density j induced 
by the applied magnetic induction B = zB. 



where fc„ = (n -h l/2)27r/W^. 

The current density j{x,y) now can be obtained from 
Eq. (3). Since V • j = 0, we also can write j = \7 x S, 
where S = zS, and S{x,y) = {B /2iJLoX^)s(x,y) is the 
stream function given by 

( \ 2,8 ^ (-l)"cosh[A;„(a;-L/4)]cos(fc„y) 
s{x,y) = y 



n=0 



fc^ cosh(fc„Z//4) 



(10) 

The series converges rapidly, and Fig. 2 shows the result 
for X > obtained by summing over n from to 10 when 

L/W = 2. 

Since Eq. (3) yields 3y{Q,y) = -((/)o/27r/xoA2)7y (0, y), 
Eq. (6) can be integrated to obtain (assuming the addi- 
tive constant A70 = 0) 

^7(2/) = ^ E tanh(fc„L/4) sin(fc„y). (11) 



ra=0 



U3 



(9) 



As shown in the Appendix, when L — >■ 00, the sum can be 
expressed in terms of the Lerch transcendent $(z,s,a), 
and in this limit A'y{y) also can be expressed as Aj{y) = 
{4:BW'^/(j)o)ipo{TTy/W), where ipoilJ^) is a function defined 
in Ref. 7. This is the explanation of why (po(/x) was found 
to be a material-independent universal function'^ when 
A » W and £ » VF. 

For all ratios of L/W, the maximum value of Aj{y) 
occurs at y = W/2, where 

14:C(3)BW^ BW"^ 

^,(wft) = --S^=,.,n?^.L«w. (13) 
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FIG. 3: A7(W/2) (solid) vs LjW and the approximations of y\G. 5: /c(B)//c(0) vs BW'^l4)o calculated from Eqs. (11) 
Eqs. (12) (dot-dashed), (13) (dotted), and (14) (dashed). and (15) for LjW = oo (solid), 1 (dashed) and 1/2 (dotted). 
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FIG. 4: A7(i/)/A7(Vy/2) vs y/(W"/2) for L ^- oo (upper 
solid curve), LjW = 1 (dashed), LjW = 0.5 (dot-dashed), 
I//W = 0.1 (dotted), and i <C W (lower solid curve), for 
which A7(t/)/A7(VK/2) = t//(VK/2). For comparison, the top 
long-dashed curve shows sin(7r?//VF). 



The solid curve in Fig. 3 shows A7(W/2) normalized to 
BW"^ jc^Q as a function of LjW along with the limiting 
behaviors of Eqs. (12) (dot-dashed) and (13) (dotted). 
The dashed curve shows the interpolating function, 



A7(M^/2) 



14C(3)BW^2 



tanh 



28C(3)W" 



(14) 



where C(3) = 1.20206 is the Riemann zeta function. 

The plots of Aj(y)/Aj{W/2) vs y/{W/2) in Fig. 4 
show how the gauge-invariant phase difference depends 
upon the ratio L/W. For L/W — > oo, the curve lies 
below sm{TTy /W), shown as the long-dashed curve, and 
for L<^W, A-f{y)/Aj{W/2) = y/{W/2), a straight line. 

The maximum Josephson current Ic{B), the maximum 
integral of jcdsmAj{y) over y from —W/2 to W/2, oc- 
curs when A70 = ±7r/2, such that 



Ic{B) 



2 

W 



w/2 



cos[A'y {y)]dy 



(15) 



where Aj{y) is given in Eq. (11). Figure 5 shows plots of 
Ic{B)/I^{0) vs for L/W = 00, 1, and 1/2. The 

stretching out of the pattern along the horizontal axis as 
L decreases is easily understood with the help of Fig. 3. 

Let us define ABi as the value of B at which Jc(-B) 
has its first zero, AB2 as the difFcrence of the values at 
which Ic{B) has its second and first zeros, and ABn as 
the difference of the values at which Ic{B) has its nth 
and (n — l)th zeros. For all finite values of L/W, the 
ABn are smaller for small n than for large n. However, 
for large n, the ABn approach the limiting value 



n=0 



tanh[(2n+ l)7rL/4VF] 
{2n + iy ■ 



(16) 



To illustrate this, if we approximate A7(y) for i — > 00 by 
A7(VF/2) sin(7r2//Vt^) (see the long-dashed curve in Fig. 
4), as in Ref. 7, then the integral in Eq. 15 can be eval- 
uated in terms of the Bessel function Jq with the result 



■M 



14C(3)gW^ 



(17) 



For L/W 00, the sum in Eq. (16) can be evaluated as 
AB = [TT^/UC{3)](po/W^ = 1.8420o/W^^ (18) 

as pointed out in Ref. 15. Using Eq. (17) and the well- 
known zeros of Jo{x), we find the following values for n 
= 1, 2, 3, 4, and 5: ABn/AB = 0.7655, 0.9916, 0.9975, 
0.9988, and 0.9993. However, to evaluate the zeros of 
Ic(B)/ lc{0) without using the Bessel-function approxi- 
mation, we must numerically evaluate Eq. (15). This 
yields the following more accurate values for n = 1, 2, 
3, 4, and 5: ABn/AB = 0.8173, 0.9866, 0.9946, 0.9968, 
and 0.9979. 

When di L <^ W, the gauge-invariant phase differ- 
ence A7(y) of Eq. (11) becomes linear in y, and in this 
case we have the familiar Fraunhofer-like pattern, 



Je(g) _ sin(7rgW^L/20o) 
Jc(0) ~ '!tBWL/2(I)o 



(19) 
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FIG. 6: AS vs L/W calculated from Eqs. (16) (solid), (18) 
(dot-dashed), (20) (dotted), and (21) (dashed). 



IV. SUMMARY 

In this paper I have considered a Josephson junction 
bisecting a rectangular superconducting thin film of large 
Pearl length A = 2}? /d subjected to a perpendicular 
magnetic induction B. I calculated the gauge-invariant 
phase difference and used it to determine the B depen- 
dence of the Josephson critical current density Ic{B). 
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such that all the AB„ are the same and equal to 
AB = 2<l)o/WL. 



(20) 



The magnitude of this AB agrees with that in sandwich- 
typo Josephson junctions with thickness X along the 
z direction only in the limit di L < A.^^""'^^ 

The solid curve in Fig. 6 shows AB, the large-n limit of 
ABn, calculated via Eq. (16) as a function of L/W, along 
with the expressions for AB in the limits L/W oo 
(dot-dashed) and L <^ W (dotted). The dashed curve 
shows the approximate interpolating function obtained 
from Eq. (14), 



AB 



TT^(t>Q 



(l4C(3)W^2 



/ tanh 



'L 



28C(3)VF 



(21) 



Appendix A: The limit L —¥ oo 

In the limit L — >■ oo, the gauge- invariant phase dif- 
ference A'y{y) given in Eq. (11) can be expressed as 
Ajiy) = {WBW'^/iT'^^oMiTy/W), where 

OO 

a(V>) = ^(-l)"sin[(2n-l-l)V]/(2n+l)3 (Al) 

n=0 

= (i/16)e-^'^[$(-e-2''/',3,l/2) 

-e2''^$(-e=^'*,3,l/2)], (A2) 

and $(z, s, a) = YL^=o^''l{^ + ^^Y the Lerch 
transcendent. Note that cr(0) = and cr(7r/2) = 
7C(3)/8 = 1.0518. 
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